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The specific optical nonlinearities inherent in aggregates of metal nanoparticles under pico- and 
nanosecond pulsed laser irradiation are studied in nanoparticle aggregates formed in silver hydrosols. 
The results of experimental studies of the correlation between the degree of aggregation of silver 
hydrosols and their nonlinear refraction index (ri2) at the wavelengths 0.532 and 1.064 /jm are 
discussed. The experiments revealed that ri2 changes its sign at 1.064 /im as the degree of the 
hydrosol aggregation grows. The role of various processes occurring in resonant domains of aggre- 
gates and the kinetics of these processes under laser irradiation resulting in dynamic variation of the 
polarizability of aggregates are analyzed. The areas under study included the kinetics of particles 
displacement considering dissipative forces, heating of the particles and of the surrounding medium 
depending on the wavelength, intensity and duration of laser pulses. A theory of interaction of 
laser radiation with an elementary type domain - two bound silver nanoparticles - was developed to 
describe the kinetics of resonant domains photomodification in aggregates. This theory takes into 
account thermal, elastic, electrostatic and light induced effects. The experimental results on laser 
photomodification of silver particle aggregates are discussed in the context of our model. These 
results include photochromic and nonlinear optical effects, in particular, the nonlinear refraction 
and nonlinear absorption. 



PACS numbers: 42.65.-k, 78.67.Bf 

I. INTRODUCTION 

Nonlinear optical media consisting of aggregates of 
metal nanoparticles have attracted significant attention 
over the past decade due to their unique properties and 
high application potential (see, in particular— ). 

Primary examples of such media are sols containing 
fractal aggregates of colloid silver nanospheres and ran- 
dom metal-dielectric nanocomposites 1 ' 2,3 . Resonant ex- 
citation of surface plasmons in such systems results in 
a giant enhancement of local electromagnetic fields 1 -^. 
This enhancement was observed in the spectral range 
where the linear absorption spectra of aggregates were 
inhomogeneously broadened while the quality factor of 
the surface plasmon modes was found to increase with 
the wavelength. 

New physical models proposed in 1980-90's (see, 
e.gj^&L^), predicting various nontrivial optical proper- 
ties of these materials gave an impetus to further inves- 
tigations. Already the early experiments with silver hy- 
drosols containing fractal aggregates proved that these 
materials feature unique nonlinear optical properties^ 
(also see the review in 1 -! 2 -! 3 -. 

It was found in^ that the intensity of the four-wave 
degenerate mixing is enhanced by five to six orders 
of magnitude in aggregated silver hydrosols as com- 
pared to nonaggregated hydrosols. Other nonlinear 
optical effects experimentally discovered in aggregates 
of silver and gold nanoparticles include nonlinear op- 
tical activity (nonlinear gyrotropy) 10 , spontaneous ro- 



tation of the polarization ellipse 10 , nonlinear absorp- 
tion and refractio n 10 ! 11 ! 12 ! 13 ! 14 , the reverse Faraday ef- 
feel—————, optical Kerr effec t 12 ' 15 , harmonic gener- 
ation 16 , enhancement of nonlinear responses of organic 
molecules adsorbed on nanoparticle surfaces^ 7 - and pho- 
tochrome reactions that mediate the effect of optical 
memory 1 ^! 20 ! 21 ! 22 ! 23 . 

A considerable wealth of knowledge about nonlinear 
optical effects in random heterogeneous nanostructures 
has been accumulated by now. However, there is still a 
lack of a comprehensive understanding of the origin and 
the mechanisms of various nonlinear optical effects, espe- 
cially for fast nonlinearities with the characteristic times 
of the order of 10~ 12 s. It should be noted that non- 
linearities with the times from 10 -14 s to 10 -12 s due 
to optical excitation of electrons in individual nanoparti- 
cles were discussed in Refsi 24 ' 25 . As for the mechanisms 
of slow nonlinearity (with the characteristic times of the 
order of 10~ 7 s), these are usually associated with ther- 
mal expansion and electrostriction and are outside the 
scope of our discussion. 

One of the sources of fast nonlinearity is the effect of 
spectrum- and polarization-selective transparency in ag- 
gregated sols induced by strong laser pulses. This ef- 
fect was discovered in Rcf.— and studied in detail in 
Refsi£ ' 19 i 20 ' 22 . It happens when the laser frequency falls 
within the spectral range of resonant surface plasmon 
excitation of aggregates and the laser power exceeds a 
certain threshold. The effect is attributed to the local 
photomodification of aggregates and has the characteris- 
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tic times from 60 ps to 150 ps^i i 9 i 20 i 21 depending on the ir- 
radiation intensity. Photomodification of the nanoaggre- 
gate structure takes place at laser intensities going above 
a certain threshold. However, the photomodification is 
confined to relatively small resonance domain a 26 i 27 ' 28 and 
does not affect the large scale geometry of the aggregate. 

The optical memory effect is associated with laser pho- 
tomodification of nanoaggregates. It is believed that the 
structure of resonant domains of aggregates is subject 
to selective changes during this process^. The changes 
take place both in the size, shape and state of the reso- 
nant particles (their fusion or evaporation) as well as in 
their relative arrangements. Thus, for example, in Ref. 15 
it is suggested that we deal here with dynamic plas- 
mon resonances of silver nanoparticles whose frequency 
changes under the applied pulsed radiation. According 
toi^, this is caused by evaporation of particles in the 
laser pulsed field and the formation of a rapidly expand- 
ing plasma inside which the electron concentration is de- 
creased. I n 10 i 29 i 30 ' 31 i 32 a possibility of particles displace- 
ment in aggregates is shown under the action of light in- 
duced dipole (generally - multipole) forces. The collapsed 
adsorption layer of particles during photo- and thermo- 
emission of electrons from the metal core was shown to 
result in the same effect as in photostimulated aggrega- 
tion of metal sols^ 3 -. 

In order to analyze the regularities revealed in non- 
linear refraction and nonlinear absorption of metal hy- 
drosols, one should consider in detail the whole set of 
processes of interaction of aggregate particles with the 
external radiation, including the processes occurring at 
the interface. These regularities include, in particular, 
the change of the nonlinear refraction sign with the grow- 
ing degree of aggregation of silver hydroso l 15 ' 16 at 1.064 
fjan. However at present the physics of photomodifica- 
tion of fractal nanoparticle aggregates with plasmon ab- 
sorption has not been studied in enough detail yet and 
a theory that could be used to describe the kinetics of 
photomodification and to determine the role of the pro- 
cesses underlying this effect is not available either. Such 
a theory would help to unriddle some of the perplexing 
issues concerning the nature of photomodification and 
the nonlinear optical phenomena observed in sols with 
fractal aggregates under applied laser fields. 

The aim of this paper is to study the mechanisms of 
fast (with times of order 10~ 8 — 10~ 12 s) nonlinear-optical 
responses of nanoparticle aggregates in metal sols (above 
all, in the sols with a liquid surrounding medium) induced 
by the optical Kerr effect which can help to explain the 
features of nonlinear optical processes in metal sols. 



II. EXPERIMENTAL RESULTS AND 
DISCUSSION 

Silver hydrosols are the most suitable model medium 
for studying optical phenomena in metal colloids. Their 
optical properties in the visible range are determined by 



an isolated surface plasmon resonance lying outside the 
interband absorbtion range. Coagulation of particles in 
the process of their Brownian motion leads to formation 
of colloid structures with a fractal distribution of parti- 
cles (see Refs. ini^i 3 .). 

The growth of fractal aggregates results in a strong in- 
homogeneous broadening of the plasmon absorption spec- 
trum of hydrosols due to a strong electrodynamic inter- 
action between the particles^ 3 -. Quite a number of ex- 
perimental and theoretical papers^ 3 -^ have proven by 
now the existence of a correlation between the degree of 
aggregation of colloids and their optical properties. 

The simplest preparation technique of a standard Ag 
hydrosol uses collargol (a silver antiseptic), containing 
disperse silver (the methods of preparation of sols are de- 
scribed, in particular, ir* 9 ^). The sizes of silver nanoag- 
gregates in this type of sols are within the range 100- 
1000 nm, the average size of silver particles being 4 to 20 
nm and a volume fraction of silver being of the order of 
(1^5) • 10~ 6 . 

Particles of this type of hydrosol feature a water solu- 
ble polymeric adsorption layers (adlayer) which prevents 
their fast coagulation and direct contact of particles in 
aggregates, as well as enhances the sedimentation stabil- 
ity. 

In the used type of hydrosol the adlayers consisted of 
casein, the adlayer thickness being 10-20% of the par- 
ticle radius. In the absence of an adlayer, a hydrosol 
with the given concentration of the disperse phase would 
aggregate within less than 10~ 2 s followed by fast sedi- 
mentation. The thickness of the adlayer particles appears 
to be several times smaller in aggregates than in isolated 
particle o 3 ' 33 . In the case of collargol, the thickness of the 
adlayer was reduced to accelerate the aggregation process 
by adding NaOH (0.1 mol/1) into the solution, which 
leads to dehydration of the casein shell and reduces its 
volume. 

The degree of aggregation (A) was monitored by the 
absorption spectra of hydrosols. Measuring of the ab- 
sorption spectrum broadening to determine the degree 
of aggregation was suggested, in particular, ir* 3 -^ (where 
A is the normalized broadening of the spectrum in the 
range ^ A ^ 1). Actually, the proposed parameter, 
A, is defined as the difference between specific integrals 
over the functions describing one of the spectral curves 
at the stage of aggregation (2 — 4 in Fig. [T|) and the 
original Curve 1 in the wavelength range from approxi- 
mately Ai > 440 nm (the first point of intersection of the 
curves) to the long-wavelength point of intersection A2 

(A oc^ Ef s (X)dX ~ j£ Ef s (X)d\) 1 Ef bs (X) is the spec- 
tral form-factor (contour) of the i-th absorption band. In 
particular, the A value corresponding to Curve 3 in Fig.[T] 
is equal to 0.8. In&SL it was shown that the parameter 
A correlates with the rate of coagulation, i.e. A is pro- 
portionate to Ao(£) -1 (where No(t) is the full number 
of particles in the system including many-particle ones) 
over the time of double reduction of the initial (N) con- 
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FIG. 1: Typical absorption spectra of collargol-based silver 
hydrosols for various degrees of aggregation A: Curve 1 is a 
nonaggregated sol, Curves 2-4 show the hydrosol at different 
stages of aggregation. Curve 4 corresponds to the developed 
stage of hydrosol aggregation. The absorption spectra of silver 
hydrosols differ from each other at various stages of aggrega- 
tion, in particular: 1 - A—0 (isolated particles), 2 - A=0.2, 
3 - A— 0.8 (intermediate stages of aggregation), 3 - A—l (all 
particles are incorporated in fractal aggregates). 



centration N (t) = N/2. 

In most cases, single beam longitudinal scanning 
(along the designated Z axis)^ is a convenient method 
to study nonlinear optical proiperties of materials. The 
advantage of this method is in its simplicity and high sen- 
sitivity of measuring nonlinear refraction and absorption 
indices. We used a Q-switched Nd : Y AG laser in our 
experiments. The behaviour of the nonlinear refractive 
index ni was studied at the laser wavelength A= 1.064 
^m (W^=12mJ, r=15 ns). 

The experimental setup was built using the arrange- 
ment with a limiting apertur o 15 ' 16 ! 36 . A 1 mm diameter 
aperture, which transmitted ~1% of laser radiation, was 
placed at a distance of 100 cm from the focal area. Be- 
hind the aperture, a photodiode was placed whose output 
signal was measured. To avoid the influence of instabil- 
ity of the laser output on the measurement result, the 
signal detected with the photodiode was normalized to 
the indications from another diode placed in front of the 
focusing lens. The laser radiation was focused with an 
18-cm focal length lens. 

A 2-mm thick quartz cell containing silver hydrosol 
travelled along the optical axis Z driven by a micropo- 
sitioner, passing the areas in front and behind the focus 
through the focus. The focal spot size at the beam waist 
was 100 /Ltm and the highest radiation intensity achieved 
was 8xl0 9 W/cm 2 . The arrangement with a limiting 
aperture enabled us to measure both the sign and the 
absolute value of the nonlinear refractive index ri2 of the 
medium under study. 

Fig. [5] shows typical dependencies T(Z) of the Ag hy- 
drosol cell transmission on the cell position with respect 



FIG. 2: Normalized transmission vs Ag hydrosol cell position 
in the limiting diaphragm arrangement. Curve 1 was obtained 
at the degree of aggregation ^4=0, Curve 2 - A—l. 
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FIG. 3: Nonlinear refractive index vs the degree of aggrega- 
tion. 



to the laser radiation focus. These dependencies were ob- 
tained for the initial (Curve 1: A=0, the intensity in the 
focal area is I ~ 10 9 W/cm 2 ) and the final (A=l, I ~ 10 8 
W/cm 2 ) stages of aggregation. It has been found that the 
sign of n 2 changes with the increasing A^. Similar re- 
sults were obtained in aggregated silver hydrosols when 
irradiated by single picosecond pulses t=35-50 ps 16 . 

Nonlinear optical parameters of the medium were de- 
termined from optical transmission using the Z-scan 
equatior 
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where AT p _ v is the normalized difference between the 
maximum and minimum transmission in the dependence 
T(Z), Auq is the nonlinear addition to the refractive in- 
dex (An) at the focus (Z = 0), S is the geometrical 
transmission of the aperture (a fraction of radiation in a 
diverging Gaussian beam incident on the photodetector 
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through the aperture in the absence of nonlinear effects), 
A is the radiation wavelength; L is the sample length, a 
is the linear absorption coefficient. 

For cubic nonlincarity, the refraction index n is re- 
lated to ri2 through the known equation-^* 3 ^ n = no + 
0.5n2|i?| 2 = no + An, where no is linear refractive index, 
E is the amplitude of field intensity. The nonlinear re- 
fractive index of silver hydrosol calculated from Equation 
© 

An important experimental result obtained in this pa- 
per, which should be discussed separately, is the change 
of the sign of with increasing A (Fig. [3]). In fact, the 
initial ability of the medium for self- focusing transforms 
into its ability for self-defocusing with the increasing A. 
Explanation of this effect gets our primary attention in 
this paper because it gives us an instrument to estimate 
the accuracy of our model and whether this model can be 
used for explanation of other nonlinear phenomena. The 
first attempt to explain the effect was made ir>i^; some 
of the important questions however (see below) were left 
unanswered. In our paper as well as ir>i£ it is assumed 
that nonlinear refraction in Ag hydrosols is caused by the 
optical Kerr effect at resonances existing in the medium, 
because, as follows from our estimates, the contribution 
from thermal expansion of the medium is negligible in 
the case of nanosecond pulses used in the experiments. 

The nonlinear refractive index of a medium is re- 
lated to the third-order nonlinear susceptibility as ?i2 oc 
Uq 2 Re[x^^- The Kerr nonlinear susceptibility (its real 
part) i?e[%( 3 )] of a multicomponent system for the pro- 
cess lu = uj + uj — lu in the case of a simplest two-level 
system away from resonances can be described as (see, 
for example, (see, for example^) 

Re[ X (3 \uj,~uj,uj)} ex Y.NiidutflMi ~ cp 3 , (2) 

where (di2)i is the matrix element of the electric dipole 
moment at the frequency (^12)4, iVj is the concentration 
of particles of the ith component of the medium (sum- 
mation is performed over all components). One can see 
from this expression that the sign of ni is determined by 
the sign of the resonance detuning. According to exper- 
imental data the nonlinear refractive index of water is 
close to zero. Nevertheless it should be noted that non- 
linear dispersion properties of water can change under 
the influence of nanoparticlcs. 

In out paper as well as ini£ we assume that nonlinear 
refraction in Ag hydrosols is caused by the optical Kerr 
effect at the resonances existing in the medium, because, 
as follows from our estimates, the contribution from ther- 
mal expansion of the medium in the case of nanosecond 
pulses used in experiments is negligible. 

The fact of changing of 712 for hydrosols upon their ag- 
gregation is known; however, the models describing this 
effect do not predict the sign change (at least, upon ir- 
radiation at 1.064 /mi). Therefore, this effect should be 
additionally analyzed. 

According to experimental data nonlinear refractive in- 
dex of a water at 1.064 fj,m is close to zero. Nevertheless 



we note that nonlinear dispersion properties of water can 
change under the influence of nanoparticles. 

The primary idea of this paper is based on the assump- 
tion that the local structure of nanoaggregates in the res- 
onant domain area changes under the applied laser field 
due to the changes in the distances between adjacent par- 
ticles. We shall call this change as " photomodification" 
in our further discussion. The change of the local struc- 
ture causes changes in electrodynamic interparticle inter- 
actions, which affects the plasmon absorption spectra to 
result in the change of the denominator sign in Expres- 
sion {2J . Such local reconstruction can be associated with 
the shift of particles both under the influence of light in- 
duced multipole interaction and due to the collapse of 
adlayers of the particles, as shown in 3 - 3 -, and it is believed 
to be one of the reasons for photostimulated aggregation 
of metal sols. What would be the consequences of irra- 
diation of nanoaggregates by high intensity laser light? 



A. Optical Kerr effect in the case when 
photomodification of aggregates is negligible 

Employing Expression © for qualitative description, 
one can see that the sign of n-i depends on the resonant 
detuning sign. If the condition uj < uiyi is fulfilled, self- 
focusing is observed in a one-component medium, while 
for lo > lj\2 self-defocusing is observed. If upon summa- 
tion over components of the medium the partial contri- 
bution of some of the components dominates, the sign of 
ni will also depend on {d\2)i and the concentration of the 
ith component. It is quite possible that the dominating 
role of one of the resonances in the alteration of nonlin- 
ear dispersion properties of the medium and interference 
of the Kerr nonlinear polarizations at certain resonances 
are the factors resulting in a positive or negative nonlin- 
ear refraction. In our case, silver nanoparticles and water 
were the main components of the medium. 

Consider contribution from each of the main compo- 
nents of hydrosol into i?e[x < - 3 - ) ] at the laser wavelength 
1.064 fxva taking into account the resonances typical for 
such kind of media. In the case of silver particles, a 
surface plasmon with the absorption band in the range 
A p z =0.4-0.42 fj,m and the interband absorption band ly- 
ing in the range A < 0.35 /im play the main role in the 
formation of a nonlinear response. The frequencies of 
both these bands are higher than the laser frequency (toi) 
hence their contribution to nonlinear refraction will be 
positive. The inhomogencous broadening of the surface 
plasmon spectrum along with the effective red shift of the 
absorption band (Fig. [JJ during aggregation of silver par- 
ticles results in the increase of the positive contribution 
to nonlinear refraction. 

However just the aggregation of Ag particles on its own 
should not result in the change of the nonlinear refrac- 
tion sign because even in the case of maximum spectral 
broadening and partial absorption in the region uj < u>i, 
the weight contribution to of the resonances located 
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in the region x^ 3 ) prevails. 

The Kerr susceptibility i?e[x^ 3 '] of water at A = 
1.064 [Mm is associated with the influence of the clcctron- 
vibronic band A < 0.19 /im) having the integrated 
{dm)e-v 10~ 18 esu and of the vibronic band having the 
fundamental transition at (A «2.9 /im) and its first over- 
tone at (A » 1.45 fim)2&22.. 

The contributions from these bands to x^ have oppo- 
site signs, and the net nonlinear refraction will be deter- 
mined by the {d\2)e-v to [d\2) v ration. Our estimates 
show that the contribution from the short-wavelength 
transitions dominates, and nonlinear refraction of water 
at A = 1.064 /xm is positive, which agrees with the ex- 
perimental dat a 40 ! 41 . 

This is due to (di2) 4 _JAw e -tj| -3 somewhat higher 
than (c?i2) 4 | Aw^l -3 (Aw e _„ and Aoj v are the detunings 
of the respective resonances from u>i). 

That is, the contribution from the both medium com- 
ponents to -Refx*- 3 )] is positive and should increase during 
aggregation of particles. This however contradicts the ex- 
perimental data and requires the model to be adjusted. 

The thermal mechanism underlying vibrational exci- 
tation of water molecules, which are in contact with 
strongly heated silver particles or with the resultant 
plasma, is another probable cause of changes occurring 
in the hydrosol nonlinear dispersive properties. This 
leads to the population of high vibrational states of water 
molecules (v ^> 1), the electric dipole moment between 
which, as known, increases as {d\2) v oc y/v^ (where v is 
the number of the vibrational level). A similar mecha- 
nism of self-focusing due to vibrational rotational, and 
electronic excitation of molecules and atoms was dis- 
cussed in Ref^ 3 -. 

The estimated minimal value of unperturbed value 
{d\2)pi based on spectroscopic data for Ag sol with the 
particle size ~10 nm is 2.5xl0 -16 esu. To estimate (di2)v 
for the vibrational transition of water (v=0-l) nearest to 
uji one can use the values of the optical constants of wa- 
ter at this frequenc y 38 ' 39 and the Lorenz-Lorentz formula. 
The calculation yields (dia)^ ~10~ 20 esu. 

With the above values of (d±2) P i and (di 2 )v and the 
ratio | Aw 3 ; /A(jj 3 |ps 10 4 , and for the volume fraction of 

silver (Na 9 ^ 10~ 6 ), the ratio of nonlinear polarizations 
at the transitions in silver and in water for nonaggregated 
sols, 

p[ A 9 ■ P h]o « [N Ag (d 12 f pl \Au; pl \- 3 } : [(d 12 ) 4 |A^|- 3 ] 

(3) 

will be equal to 10 7 . 

The experimental data show that the value of x f° r 
hydrosol is several orders of magnitude higher as com- 
pared to pure water, although it is underestimated tak- 
ing into account the large value of (d\2)pi- This suggests 
that silver plays an important interceding role in the en- 
hancement of optical nonlinearity of hydrosol. 

Consider some examples describing the conditions for 
manifestation of this mechanism. The dissociation en- 
ergy of a water molecule is close to the energy of the 



upper vibrational state and is 8.4xl0 -20 J (see Ref4^), 
which corresponds to the temperature 6x 10 4 K. The en- 
ergy ~10 _10 J absorbed by a particle is sufficient for ther- 
mal dissociation of 10 8 molecules, which occupy the vol- 
ume 3xl0 -15 cm 3 and can be arranged around the par- 
ticle in a spherical layer of thickness i?=10~ 5 cm. 

The time during which all the molecules in the layer 
will be dissociated is t = R/V s ~ 10~ 10 s, where V s is 
the speed of sound in water. 

This means that the volume occupied by thermally ex- 
cited water molecules is at least 10 3 times larger than 
the volume of solver particles, although the amount Nt 
of excited water molecules does not exceed 10~ 3 — 10~ 2 . 

Assuming that (di 2 )v increases by at least one or- 
der of magnitude, and taking into account the ration 
|Acj e _„| 3 /|Ao;^| 3 pa 10 6 and the above range of values of 
Nt, the contribution of vibrational transitions of water 
into the Kerr susceptibility x*- 3 ^ can become equal to the 
contribution of the vibronic band, substantially changing 
thereby Ratio ([3]). 



B. Contribution of thermal expansion of medium 

Laser photomodification of nanoaggregates as one of 
the factors accountable for the sign change of n 2 in 
nanosecond pulses should be analyzed further from the 
point of view of kinetics of this process when the exper- 
iment is carried out with picosecond pulses. Will the 
resonant domain have time to change its structure (at 
moderated intensities) or to evaporate (at high intensi- 
ties) within the pulse duration t—30 ps? 

To answer this question, we will use the time depen- 
dencies (t p ) describing the development of photomodifi- 
cation in such Ag hydrosol o 9 ' 21 . It follows from 21 - that t p 
depends on the pulse energy density and decreases from 
150 to 50 ps upon the fourfold increase of the energy 
density. Obviously, at higher energy densities, the time 
t p can further decrease, and the particles can evaporate 
during one pulse, i.e. the role of this process in the Kerr 
nonlinearity seems to be realistic from the point of view 
of kinetics for picosecond pulses. 

Note, that the use of picosecond pulses in Z-scan tech- 
nique led to similar: the change of a sign n 2 with the 
subsequent increase of its absolute value in the process 
of aggregation of sols 16 . 

This lead us to important conclusion about the negligi- 
ble contribution of the thermal component of nonlinearity 
caused by heating of the medium under applied radiation 
and a subsequent expansion of the heated region. Such 
heating reduces the medium density and hence decreases 
its refractive index as is the case with self-defocusing. 
The negligible contribution of the thermal mechanism to 
sclf-dcfocusing is confirmed by an insignificant expansion 
of the heated region even compared to the focal waist 
size during irradiation by a nanosecond pulse. Indeed, 
the size of the thermal lens being formed is ro = V s x t, 
and for r =15 ns pulses it does not exceed ~ 10~ 3 cm. 
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For 30 ps pulses r will be only ~ 10 _6 <4, where db 
is the beam diameter. This is an order of magnitude 
smaller than the experimental size of the focal waist and 
at least two orders of magnitude smaller than the beam 
diameter outside the focus (1 cm away from the focus) 
where self-defocusing is maximumal (see. Fig. [2]). 

When analyzing the conditions for manifestation of 
nonlinear refraction, one should take also into account 
the scattering of laser radiation by metal hydrosols re- 
lated to the formation of vapor bubbles around nanopar- 
ticles absorbing radiation 4 - (the time of their formation 
and growth was 10 -8 s (see Ref4^)). 

However, when a cell with hydrosol is located in front 
of the focus this effect will only prevent the Kerr self- 
defocusing. This means that the Kerr nonlinearity on dy- 
namic plasmon resonances of silver and vibrational tran- 
sitions in water makes a much greater contribution than 
under real experimental conditions. 

In estimating dynamic frequency shifts one should also 
also allow for the shift of u> p i due to the change in 
the permittivity (e m ) of the environment when a va- 
por cavity is formed around the particle, because u p i ~ 
l/yT + 2e m Ref. 4e . This a correction does not exceed 

pi) vapor J pi) water ^1-12. 

Note that the above analysis did not take into account 
changes in the local structure of the nanoparticle aggre- 
gates themselves under modifying laser field. These can 
be done by modeling the processes and their kinetics oc- 
curring in resonant domains of large aggregates under 
the action of electromagnetic radiation. The following 
section is devoted to the formulation of such a model. 



III. DYNAMIC LOCAL RESTRUCTURING OF 
FRACTAL NANOPARTICLE AGGREGATES AS 
A DOMINANT SOURCE OF OPTICAL 
NONLINEARITY AND PHOTOCHROME 
EFFECTS 

A. External field model of resonant nanoaggregate 
domain 

Consider a model of interaction of a simplest aggregate 
- a bisphere consisting of two spherical nanoparticlcs of 
equal size, with the external radiation. Besides, the bi- 
sphere (dimer) or its modification - trimer (consisting of 
three particles), can be treated as a component of a more 
complex aggregate in which it plays the role of a resonant 
domain. Such domains that are resonant to the given 
wavelength are randomly located at different sites of the 
large aggregate. Of course, such a domain is not elec- 
trodynamically isolated in a fractal aggregate, but the 
knowledge of its basic properties and their change under 
the influence of surrounding particles (eigen frequencies, 
relaxation constants, etc) would help to predict the re- 
sponse of the fractal aggregate to the resonant laser radi- 
ation. To build such model, let us considered basic types 
of interactions in a bisphere (interaction of particles with 
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FIG. 4: The absorption spectrum of Ag bispheres in the lon- 
gitudinal direction of polarization direction and for various 
interparticle gap values h. The inset shows the bisphere struc- 
ture and the external field orientation (the bisphere particle 
radius is ro = 10 nm, V is the volume of particles) . Solid line 
corresponds to h=0 nm. 



each other and with the external optical radiation). 



1. Interaction with radiation 

We will limit our consideration to the discussion of 
the interaction of optical radiation with the elementary 
model of a resonant bisphere consisting of two identi- 
cal spherical silver nanoparticles (Fig. [TJ. The dipole 
moment of such a bisphere P, generated by field E, is 
P = ayE, the polarisability of bisphere is ay is repre- 
sented by the sum of responses of single oscillators: 

where LOoj is resonant frequency of jth oscillator, Tj is the 
relaxation constant, Aj is coefficient. And u>oji Tj, Aj 
are the functions of a particle radius (ro) and distances 
between surfaces of the particles included in bisphere (h). 

An absorption spectrum in the longitudinal direction 
of polarization has been calculated for Ag bisphere par- 
ticles with ro = 10 nm using the multipole approach 34 
(with the account of up to the 32th order multipoles). 
The spectrum has been calculated for different values 
of the interparticle gap, h, and for the medium refrac- 
tive index in the wavelength region under study being 
n=1.4 (Fig. 2]). The optical constants for silver were 
taken iromM- and adjusted for the finite size-effect in 
small metallic sphere o 46 i 48 . 

As seen in the Fig. [4] the plasmon absorption spec- 
trum (A > 400 nm) is represented by oscillating func- 
tion. To find the analytical form of this function, which 
is necessary in further calculations, it is convenient to 
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use the approximated expressions and to describe spec- 
tral curve as the sum of quazi-lorentzian contours of two 
single-resonant oscillators. The approximated analyti- 
cal expressions for resonant frequencies of the long- wave 
and the short-wave oscillators (neglecting low-amplitude 
maxima in the wing of the bisphcre absorption spectrum) 
are 



ujqi = fioi • 10 1{, s-\n i = 4.93 exp 



lo 02 = ^02 ■ 10 15 s- x ,O 2 = 4.93 exp 



0.252 
0.32 + h 



0.345 
0.72 + h 



(5) 



(6) 

where h is interparticle gap (in nm). 

As particles move apart (h — > oo) Woi approaches 
the plasmon resonance frequency uo p i of an isolated par- 
ticle (A ~ 0.39/im) whereas once W02 reaches uj p i it stops 
changing and approximation by Expression ([B]) is no 
longer valid. However considering that the bispere inter- 
action with the long- wavelength radiation (A > 0.6 /im) 
is more important in our model, the absorption band 
with wo2 in this case can be ignored. 

The values of the amplitude factors A\ , A% in were 
approximated by expressions: 



A x = 



10 4 c 



Aire 



1/2 



VTioi, A 2 



10 4 c 



Aire 



1/2 



VT 2 a' 2 , (7) 



1.28- 10 2 



1,93+ (fioi ~4.96) ; 



■(cm" 1 ), a' 2 = 280(cm" 1 ), 

(8) 

where a[ and a' 2 are absorption cross-sections of radiation 
in the unit volume of a particle at exact resonance (see 
Fig-H]), V is a particle volume, e m is the relative dielectric 
permittivity of medium (in the case of hydrosols it is a 
water), c is the speed of light in vacuum. 

Relaxation constants T\ and T 2 are determined by the 

46 

expression u : 



Ti 



r 6 + 



V F 



(9) 



where r& = 3 ■ 10 s is the relaxation constant of free 
electrons in a metal bulk (in our case Ag), Vp is the Fermi 
velocity of electrons, tq is a particle radius. 

Knowing the bisphere polarisability a d we can deter- 
mine the absorption cross-section for laser radiation at 
the frequency: 



o d = 47T— eH 2 Im[a d {oJi)}- 
c 



(10) 



Hence, the power absorbed by the bispere is , where 
is the intensity of laser radiation. Besides, the parti- 
cles get polarized in the laser field and the light-induced 



dipoles interact with each other. Hence, the power ab- 
sorbed by the bispere is W d — cr d I, where I is the in- 
tensity of laser radiation. Besides that in the field of 
laser radiation particles are polarized and light-induced 
dipoles interact with each other. The above mentioned 
parameters are sufficient for calculation of the real part 
of the bisphere polarisability. 

In the case when the angle between the field intensity 
vector E and the bisphere axis (6) is or n (see Fig. 01, 
the interaction energy and the force F em between parti- 
cles are defined by the expressio n 10 i 29 ' 30 : 

Uem = -—I ■ Re[(a d (u;, h) - 2a (w)], F em = - d J™ , 
c an 

(11) 

where ao(u)) is the polarisability of an isolated particle 
(the second term in the expression for U em , does not 
dependent on h, and it describes the energy of individual 
particles polarized in the field). Note, that under such 
field orientation (at 8 = 0) the moment of force acting 
on a bisphere is zero. 

Thus, the accuracy of our calculations and that of the 
calculations in 1 ^ is defined by the adequacy of the Claro 

mode l2^0 1 31. 

Dispersion of the surrounding medium in the spectral 
range of interest (at A=540 and 1080 nm) was taken 
into account when calculating polarisability of the Ag 
bisphere (for water the difference in dielectric permeabil- 
ity at these wavelengths is negligible and does not fall 
outside the limits 10" 2 ). 

Thus, the action of laser radiation is actually reduced 
to two important factors: the heating of particles and 
the appearance of the interaction force of induced dipoles 
(in general case multipoles) F em , that directly influences 
the movement of particles in a bisphere. There are also 
other approaches to the optical binding forces calculation 
operating between nanoparticles in optical fields (see, 
e.g. : 



.2,49^ 



2. Estimation of a radiation pressure 

Estimate also the role of the radiation pressure force 
(F pr ). It is important to note, that this force acts orthog- 
onally to the bispher axes, therefore it doesn't change a 
relative nanoparticles position. For the same reason it 
does not have effect on the bisphere optical properties. 

We estimate F pr using the expression: 



F v - P r S, 



(12) 



where S = nr^ is the radiated particle surface, P r = 21 /c 
is the radiation pressure at full reflection of the radia- 
tion. Acceleration gained by a bisphere is a = F p /2m = 
47rr 2 //[2c(4/3)7rr^/9] = 3I/(2cr p), p is the particle den- 
sity. Taking i = 10 8 W/cm 2 , r =10~ 8 m, p ~ 10 4 kg/m 3 
(for Ag) we obtain a = 5 x 10 7 m/s 2 . Accordingly, for 
the nanosecond pulse r = 10~ 8 sat the maximum exper- 
imental intensity this will lead to a shift of the bisphere 
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FIG. 5: Deformation of identical particles with radius ro when 
their adlayers are in a contact with each other (AL is an 
adsorption layer). 



by I = ar 2 /2 = 2.5 x lCT 9 m=2.5 rim. For t = 3 x KT 11 
s we have I = 2 x 10~ 14 m. If the bisphere is the part of 
an aggregate, the shifts will be much less. 



3. Elastic interaction 

To prevent fast aggregation of particles in hydrosols, 
a water-soluble polymer is often added to the initial so- 
lution before a disperse phase starts developing. In this 
case, as has been mentioned in the previous section, each 
new particle appears to be surrounded by an adsorp- 
tion layer (a polymeric shell) , having the thickness about 
/io/2 = 1-^2 nm and more (Fig. [5]). This shell reduces 
the Van-der-Waals attraction when particles collide, pre- 
venting coagulation (sticking) of the particles. The gap 
between adjacent particles in the aggregate is about 1-^5 
nm for radius ~ 10 nm. At such distances, the Van-der- 
Vaals attraction between particles (F v ) is already high 
enough for the particles to form aggregates. The key role 
in compensation of this attraction is played by the elas- 
ticity force (F e ) arising from deformation of polymeric 
shells of adjacent particles. To describe this force, we 
will use the results obtained for the problem of defor- 
mation of two balls (Fig. [5]), which is also known as the 
Hertz deformation/tangent problem^. 

The potential energy of such an elastic interaction in 
the case of identical spheres of radius r can be described 
by the expression—, 



U e = -(h -h) 5 ^ 



r + ho/2 



1/2 



(13) 



where E e is the effective elasticity module of polymeric 
adlayer, a e is the Poisson coefficient (the mean values 
~ 0.15). Since a e is very small, it will be neglected in 
our further calculations. Unlike the classical Hertz prob- 
lem, where deformation of two contacting spheres is con- 
sidered, in our case the spheres contain a metal kernel 
inside, which at least at the initial stage of deformation 
does not affect the process of deformation of the exter- 
nal surface. When using alternative expressions for the 
energy of elastic (steric) interactions of polymeric adlay- 
ers of particles (see, for example^) it is not possible 




FIG. 6: Diagram of particles capture in a potential well. 



to avoid an arbitrary choice of adjusting parameters for 
specific types of polymeric molecules. A detailed com- 
parative analysis of elastic (steric) interaction potentials 
of particles, obtained in terms of different models, will be 
reported elsewhere. Consider now the process of coagula- 
tion of two identical particles (before the forces become 
balanced) under the action of the Van-der-Waals force 
F v , the elasticity force F e = —dU e /dh and the viscous 
friction force Ff. Fig. [6] illustrates the total interaction 
energy of particles U(h) = U v + U e (U v is the Van-der- 
Waals energy) as a function of the distance between their 
surfaces h. 

The particles approach each other driven by the force 
F v . Before their polymeric adlayers come into contact 
this movement is counteracted by the viscous friction 
force Ff only and once they come into contact the elas- 
ticity force F e adds in as well. As a result of the particles 
approach a stable position of the particles in a bisphere 
with the gap h is achieved when the secondary potential 
well in U(h) = U v + U e (Fig. [6]) satisfies the following 
conditions: 



U ± >0, U 1 = W -A ] 



Umax U 7n i n ^ 

(14) 

where Wq is an initial kinetic energy of particle (Wo ~ 0), 
Af is the work of the viscous friction force in the way 
from the capture of particle in the potential wall before a 
collision with a barrier separating the secondary potential 
wall from the main potential minimum. 

In Fig. [7] the dependence E e (ho), satisfying to condi- 
tions (|14p is shown. The dependence was derived from 
the solution of the problem on particle movement in a 
potential field, taking into account the friction force and 
selecting appropriate values of the elasticity module. In 
the range 3 nm ^ ho ^ 0.8 nm for the particles whith 
ro = lOnm it is well approximated by the function 



E e (h ) = 1.67- 10 8 -exp 3.09 (h /2 - 1.62) 2 . (15) 
Here ho is measured in nm, and E e is done in N/m 2 . 
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FIG. 7: The effective elasticity module of adlayer vs the ad- 
layer thickness ho/ '2. 



For the purpose of our study we will take the advantage 
of the data available on exactly the kind of silver sols 
with a casein adlayer which were used in the experiments. 
It should be noted, that the elasticity module values in 
the contact area chosen in ([15]) for the material of the 
particles are 4-5 orders of magnitude higher compared to 
those for a free casein gel in Refi^ 2 -. This is explained by 
the fact that it is impossible, using in the calculations 
the reported values of E e , to account for the existence 
of nanoparticles aggregates with the distances between 
their surfaces h ~ 0.5 -j- 2 nm for which the Van-der- 
Waals forces would counterbalance the elasticity force. 

Large values of the adlayer elasticity module can be 
explained by the following factors. Adlayers of contact- 
ing particles in aggregates already are in a deformed 
state near the contact surface. Besides, the adlayer is 
attracted to the particle surface due to electrostatic, Van- 
der-Waals, ion-dipole and other types of interactions. 
This all means that the adlayer is in a compressed state, 
which increases its density in comparison with a free gel. 
Moreover, the gel structure changes as well due to inter- 
action with the surface because the energy of binding of 
the particle surface with polymeric molecules turns out to 
be larger than the binding energy between the molecules 
(the latter accounts for the elasticity of a free gel). The 
thinner the polymeric adlayer, the more it is compressed 
and the greater is the contribution of interaction with 
the surface, which should result in an increased module 
of elasticity. The bonds in a polymeric gel grid are de- 
stroyed during heating of particles and adlayers induced 
by the applied laser field, which decreases the elasticity 
module of the gel. This, in turn, alters the conditions of 
the balance of forces for bisphere particles since the elas- 
ticity force decreases. The latter leads to degradation of 
the secondary potential well and the particles start ap- 
proaching to come into almost a full contact, i.e. to get 
into the main potential minimum. Thus the temperature 
dependence of the elasticity module is an important con- 



TABLE I: The potential barrier for different relaxation times 



n, s 


60 


600 


3600 


U f , eV 


1,04 


1,12 


1,17 



sideration determining the accuracy of our model. For 
example, the elasticity module E e changes from 1.3 • 10 4 
N/m 2 to 3 • 10 3 N/m 2 when the temperature is raised by 
30 A' 52 . In our case, however, the experiments did not 
reveal this kind of a strong dependence: the structure 
of silver aggregates did not show any noticeable changes 
when heated up to 370 K, which is supported by the data 
on the absorption spectra. This means that the elastic 
characteristics of polymeric adlayers of particles in ag- 
gregates are subject to a very little change, which can be 
attributed to their strong initial compression and higher 
binding energy with the particle surface. 

The rate of heating is another important factor affect- 
ing the elasticity module of a polymer in the process of 
heating. As is known from Ref^ 3 -, destruction of a poly- 
meric grid occurs within terminal time T<j. So, if it does 
not take place within the time equal to the laser pulse du- 
ration there will be no changes in the elasticity module. 
According to the theory^, an average waiting time for 
local energy fluctuation Ej (in a single atom or a small 
group of atoms bound into a grid point) depending on 
the temperature is found as 

T f ~ to • exp(E f /kT), (16) 

where To ~ 10~ 12 -j- 10 _13 s. The bond with the potential 
barrier Uf will be destroyed when Ef ^ Ut. 
The average waiting time for this to happen is 

T b ~T Q -exp(U f /kT), (17) 

which is the time required for a polymer to relax into 
a new stable state (with a smaller number of bonds be- 
tween the molecules) corresponding to the changed tem- 
perature. 

Proceeding from the fact that fractal aggregates re- 
main stable when heated up to 100°C (~ 370 K), we 
estimate Uf. The values of Uf obtained for various re- 
laxation times at ~ 370 K and To = 10~ 12 s are sum- 
marized in the Table 1. 

Assuming the average value Uf = 1.1 eV (which, inci- 
dentally, is in a good agreement with a similar value for 
polymeric media), we find the critical temperatures (T cr ) 
below which the times t& = 10 -8 and 10 -11 s are not long 
enough for relaxation to occur. We obtain t& = 10~ 8 s — 
T cr = 1385 K; r b = lO" 11 s — T cr = 4070 K. Theoret- 
ically, once these temperatures are reached, the polymer 
elasticity module should drop dramatically. In our opin- 
ion, the change of E e with temperature can be assessed 
with more accuracy in terms of the fluctuation frequency 
v r : 

v r = r^ 1 = Tq- 1 e X p{-Uf/kT). (18) 



10 



In this case it is reasonable to draw an analogy with 
the chemical reaction at the rate v T . Then the change of 
E e can be expressed as 



a) 



dEe 

dt 



— —v r E„ 



(19) 



It is easy to check, that for a short time At 
relative change satisfies the equation: 



AE e 



At 

TO 



exp(-C///fcT). 



(20) 



If At ■ Tq 1 ■ exp[-U f /kT] < 1, then AE e < E e , i.e. 
there is virtually no change. For this purpose the condi- 
tion T < U f /\n(At / To) should be satisfied which corre- 
sponds to the definition T cr . 

Note that the properties of a polymer can be strongly 
modified by other factors as well. For example, at the 
particle surface temperature T > 600 K, explosive vapor- 
ization takes place in the adjacent layer of overheated wa- 
ter with a resultant dehydration of the polymeric layer 55 . 

It is therefore necessary to take into consideration the 
whole variety of factors affecting the elastic interaction 
of particles in a bisphere. Besides, the repulsive interac- 
tions between particles are also caused by the interaction 
of electric double layers (EDL) consisting of electrolite 
ions which are adsorbed by particles along with poly- 
meric molecules. 



4- Electrostatic interaction 




FIG. 8: Electrostatic interaction of two spherical nanopar- 
ticles: a) overlapping of double electric layers (EDL) when 
particles collide; b) the structure of EDL. 



Electrostatic interactions of spherical nanoparticles 
in a bisphere are generally described in terms of the 
Derjaguin-Landau-Verwey-Overbeek (DLVO) theory (see 
for example 56 ). According to this theory, each particle 
in electrolite solution is surrounded by an electric dou- 
ble layer consisting of electrolyte ions. Electrostatic in- 
teraction is virtually absent when the distance between 
the surfaces of two particles exceeds the total thickness 
of these two layers (2Ao), and at shorter distances these 
layers overlap (Fig. [S^,) . Due to a partial displacement of 
electric charges from the external part of EDL (a counte- 
rion layer) , there is an electrostatic repulsion of the inter- 
nal, dense parts of the EDL particles (layers of potential- 
determining ions). Since concentration of electrolite ions 
can be quite high (up to 10~ 3 -f -10 -1 mol/1) in polymer- 
containing silver hydrosols, the particle surface adsorbs 
polymers as well as ions which form EDL. The thickness 
of the dense part of EDL is defined by the ions concen- 
tration and depending on this concentration it can either 
exceed the thickness of the polymeric part of the adlayer 
or be localized inside that part. The latter corresponds 
to our experimental conditions (Ao ^ 1 nm)); electro- 
static repulsion takes place only in case of appreciable 
deformation of polymeric shells when the interparticle 
gap reaches the total thickness of the dense parts of EDL 
of neighbouring particles. 



The potential energy of Coulomb interaction U q (h) of 
overlapping EDLs of spherical particles of radius ro is 
defined as (see, for example 

U q {h) = 2ire a er p 2 s In [1 + exp(-/i/A ]) , (21) 

where So is dielectric constant, e is a relative dielectric 
permittivity of the surrounding medium, Ao is the Debye- 
Hiickel radius, tps is the potential at the border of the 
Helmholtz layer (the layer thickness doesn't exceed the 
size of hydrated counterions constituting the layer). The 
latter is the parameter describing the EDL on a particle 
surface. 

Unlike the Youkawa potential (see, for example^), 
Expression (|2"Tj) is more convenient for estimation of 
U q (h) as it contains the parameter tpg. The potential 
tps for a low electrolite concentration is close to be equal 
to electrokinetic (£) potential defined experimentally by 
means of electrophoresis 5 ^ 5 -^ (the potential at the glide 
plane of a nanoparticle moving in an applied electric 
field). A further advantage of using this expression is 
that it does not require data on the effective screened 
charge of particles in the layer of potential-determining 
ions. 

For estimation of the EDL parameters, described in 
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the Stern theory (see for example 56 ), consider its struc- 
ture. In Fig.[5j3 characteristic EDL borders with a poten- 
tial of the Helmholtz layer tps and Ao (the characteristic 
extension of the diffuse part of a counterions layer) are 
shown (at the distance where the potential of the layer of 
potential-determining ions located on a particle surface 
falls by e times). In order to estimate the EDL param- 
eters described in the terms of the Stern theory (see for 
example 56 )), consider the structure of EDL (Fig. [SJd). 
Typical EDL boundaries with the Helmholtz layer po- 
tential ipg and the typical size of the diffuse part of a 
counterions layer Ao are shown in Fig. [SJd. 

In the case of spherical particle ips can be found as 5 ^: 



<P5 



<7o 



+ 



ineoero 4ireoe(r + S) ' 



(22) 



where qo is the total charge of potential-determining ions 
adsorbed on the particle surface, q c is the charge of coun- 
terions in the Helmholtz layer. 

According to the Stern theory, it's determined by 56 : 



1 + x- 1 exp(A t /T) 



-4, 



F ■ z ■ ips + $ 
R ' 



(23) 



Here qoo is the charge limit for counterions in the 
Helmholtz layer q^ ^ qo, x is the dimensionless pa- 
rameter equals to the mole fraction of counterions in the 
solution, F is the Faraday constant, z is the counterion 
charge (z=l), <E> is the energy of specific adsorption, R 
is the absolute gas constant. 

Taking into account that 5 <C r (5 is of the order of 
the counterions sizes) and assumung ~ qo, we obtain 
from {32]) and 



'•ps 



<7o 



exp [(F-z-(p s + §)/RT] 

47re er x + exp (F ■ z ■ <p s + $) /RT] ' 



(24) 



Values qo and <& can be derived from the results of 3 -^, 
where for the same type of hydrosol at the given elec- 
trolite (NaOH) concentration (c m ) and T — 300 K 
the following values were obtained: £ = —0.064 V at 
c m ±2 10" 3 mol/1 and C = -0.034 V at c m - 10~ 2 
mol/1. Bearing in mind that ^-potential at low values of 
c m is equal to y> j 56 i 59 we can use these values to estimate 
qo and $. Using (2U, we find $ = -2.02 • 10 4 J/mol, 
qo = 9.8 • 10 -18 C, which agrees with other estimates of 
the charge of a particle 3 -^ 3 -. 

Moreover, an estimate for Ao at c m ^ lO^ 1 mol/1 and 
T = 300 K is available 3 -^, which is A ~ 10" 9 m. Since 
A oc VT, we obtain the expression for U q {h) with quite 
specific values <ps(T) and Aq(T) at c m — 10 _1 mol/1 



ips « -0.11[1 + 6exp(-1.16 • 10V/T)]-Hn f2 - 
A ~ 1.7510 • 10- 2 y ^r/c m ~(nm). { j 

In particular, at T = 300 K we obtain ips ~ —0.0108 
V for c m = 0.1 mol/1, corresponding to the experimental 
conditions in the used silver hydrosols. 



Expression ((25)) is valid for a steady EDL. In the pro- 
cess of photomodification the particles and surrounding 
medium get heated under the action of laser radiation. 
The EDL parameters change with increasing the temper- 
ature. A question arises whether EDL will have enough 
time to achieve a steady state during a heating? 

Let us estimate the time of achieving a balance of 
charges in the layer of Ao thickness: t(\q) ~ Ao/vi, v; is 
an average velocity of counterions in the changing electro- 
static field of a particle ( q is assumed to be independent 
of T). In turn 



Vi~Bi|E|, E^E S 



(qo - q c ) 

47T£ £ro 2 ' 



(26) 



where Bi is the counterions mobility in the medium, E is 
an intensity of electric field generated by the Helmholtz 
layer of ions in the diffuse part of EDL, Es is the electric 
field intensity at the external boundary of EDL. 

Taking into account (|24| . we have an estimate for Es 



\Es\ 



TO ' 



> 



Bi\ips\ 



(27) 



At r = 10~ 8 m, tp s - 0.01 V and A ~ 10~ 9 m (for 
c m = 0.1 mol/1) we obtain \E S \ ~ 10 6 V r A > 2-10 8 s for 
a mobility of counterions in water Bi ~ 4 • 10 -8 m/ (V-s). 
At the same time laser pulses used in experiments have 
durations ~ 10~ 8 and 10 -11 s. 

Thus, the structure of EDL changes very little under 
applied laser field when the temperature of the surround- 
ing medium changes. Therefore in our further discussion 
the properties of EDL particles will be assumed to re- 
main unaffected by laser radiation. Moreover, we should 
also bear in mind that the presence of a polymer in the 
adlayer can severely decrease the mobility of counteri- 
ons. Hence the time required for EDL restructuring will 
increase. 



5. The Van-der-Waals interaction 

Dispersion forces of the Van-der-Waals attraction are 
one of the key factors determining movement of particles 
in a bisphere during laser irradiation. We will employ the 
Hamaker and de Boer theory (see, for example ; 57 ' 58 ) to 
describe these forces. According to the theory, the inter- 
action energy of two spherical particles with the radius 
tq is described by the following expression: 



u v — „ 

6 



2r 2 



r 2 
'12 



4r 2 



(28) 



where r 12 = 2r + h is the distance between the centers 
of particles, Ah sa 50fcT is t ne Hamaker constant for 
metal hydrosol^ 1 -, In the calculations, a limitation is made 
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on the value here To = 300 K, k is the Boltzman con- 
stant. In calculations the limitation is made on the value 
(rf 2 - 4r 2 ,)/4r 2 l > 0.0001. Thus, the Van-der-Waals force 
acting on each particles is equal to F v = —dU v /dh. 



6. The force of viscous friction 

When particles move in a bisphere they experience 
counteraction from the medium (in hydrosols - from wa- 
ter) - the force of viscous friction. This force is estimated 
by the Stocks formula: 



F f = 6Trr]r V p , 



(29) 



where 77 is the viscosity of a medium, V p is velocity of 
particles in a bisphere relative to its center of mass. 




T 



7. Heating of particles with consideration of heat emission 
into surrounding medium 



FIG. 9: Heat exchange between the particle and the surround- 
ing medium 



The primary effect of absorption of laser radiation by 
metal particles of a biosphere is an increased internal en- 
ergy of the electron gas. Energy transfer from electrons 
to the metal crystal lattice occurs rather slowly because 
of the big difference in the mass of electrons and ions. 
Equilibrium times in the electronic r e and ionic Tj subsys- 
tems are much shorter than the equilibrium time between 
electrons and the lattice. So it is possible to introduce 
two subsystems: electronic and ionic. For r e and we 
obtain the following estimates: 



((J q 



N)" 



K-Vi-N)" 1 , 



(30) 



where a q is cross-section of Coulomb collisions, v c , vj 
are the average electron and ion velocities, respectively, 
N is the concentration of particles. 

With the typical values of of electron (ion) concentra- 
tion in a particle being N ~ 10 28 m -3 , for kT ~ 1 eV we 



have o q ~ 10 nr 
Using (|30|). we find 

r e ~ 10 



v c = 6 • 10 5 m/s, Vi — v cV /m c /mi. 



17 s; 



10 



-14 



10 



-15, 



That is equilibrium times are much shorter in com- 
parison with pulse duration and typical heating times 
Th ^ 10 -12 s (under our conditions), so it is possible to 
introduce temperatures and T e for electronic and ionic T 
subsystems, respectively. Under these conditions we can 
write down the following equations for T e and T^ 2 - 



c e — = -g{T e - T) + — , 
at v d 



where C e , Ci are electronic and lattice thermal capacities, 
respectively, g is the rate of energy exchange between 
subsystems, v d is the bisphere volume. These values are 
taken from^ 3 -, where 

C e = 68 • T e J-m^-K" 1 ; d = 2.5 • 10 6 J-m" 3 ^- 1 ; 
,g = 4-10 16 J-m^-K- 1 ^- 1 . 

The term crl/vd account for energy gain by the elec- 
tronic subsystem (to unit volume) due to absorption of 
laser radiation with intensity /, qi is heat loss from a 
particle into the environment. 

Definition of qi is rather difficult problem for which a 
non-stationary problem of heat conductivity has to be 
solved first with consideration for the phase transition 
" liquid- vapour" 55 i 64 . This essentially complicates under- 
standing of the particles kinetics in a bisphere under laser 
radiation heating. Therefore we will use the following 
simple model. The heat loss per unit volume is described 
as 



j -S aa • grad(Ti) • S 



Vd 



(32) 



where S is the surface of a particle, j is the heat flux 
density, ae is the heat conductivity factor of water. 

Let L is the boundary of heated area in surrounding 
liquid (a water) (Fig.[9|). 

The speed of movement of this border Vh ~ \J Om/t, 
/cmoPm is a thermal diffusivity coefficient of 
water, c mo , p m is a specific thermal capacity and den- 
sity, respectively. 

Then the border r of a heated area is defined by the 
equation 



Ci 



dt 



7(T e - Ti) - qi, 



(31) 



dr 



1 



^ = VT' r = r + -Va—t. 



(33) 
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Now we can estimate grad(T) : 

grad (T) « ^—^ = 2 /^L^l (Ti - T ), (34) 

where To is the temperature of water near the border of 
heated area (T = 300 K). 
Using ([32]) we find 

qi = -(ae m c mo p m ) 1/2 ^ 1/2 (T - T ). (35) 

In the model under discussion the parameters of a 
medium are assumed to be invariable, although they de- 
pend on the medium temperature, and also change dur- 
ing the phase transition "liquid-steam". However, as fol- 
lows from££, phase transition can occur at T w 600 K. 
That means that the liquid although overheated main- 
tains its aggregative state. 

8. Kinetic equations for a bisphere in the process of 
photomodification 

As follows from the above discussion, the change of the 
bisphere characteristics under the action of optical radi- 
ation, which we will call " photomodification" , depends 
on many factors such as the laser pulse energy and du- 
ration, the radiation frequency and also the initial state 
of the bisphere (the interparticle distance). The shape 
of the potential energy surface shown in Fig. [10] (for the 
laser radiation intensity I = 10 7 W/cm 2 and the initial 
distance between particle surfaces h\ = 2 nm) can give 
an idea of the movement of particles in a bisphere. 

The potential surface in absence of radiation is defined 
by the Van-der-Waals, elastic and electrostatic interac- 
tions. In Fig. 110a the surface section with to ^ 3.6 • 10 15 
s _1 describes a surface in absence of radiation (Curve 1 
in Fig. [TUb) where the potential wall is observed. The 
fold in the field cj>3.8-10 15 s _1 is the result of the light- 
induced particle interactions. These interactions, as is 
shown in Fig. [TUb . can lead to a shift of the potential 
well (Curves 2,3) or to its disappearance. That is, the 
particles which were in the stable position h± start mov- 
ing to a new stable position under the action of radiation. 
Their resonant frequency will be shifted as well both to 
the low-frequency range (Curves 2, 4 in Fig. |10b) and to 
the high-frequency range (Curve 3). 

The situation, in actual fact, is more complicated. The 
particles and the polymeric layer are heated up by laser 
radiation to result in the change of the module of elas- 
ticity and of the shape of the potential surface. The 
potential surface shown in Fig. 1101 corresponds to the be- 
ginning of influence of laser radiation and it can describe 
the movement of particles at the initial stage. In order 
to describe the bisphere behaviour (the interparticle dis- 
tance, resonant frequencies, temperature), it is necessary 
to solve the system of kinetics equations in which the 
factors analyzed in the previous sections are taken into 
account. This system can be written down as follows: 



— — 2V nV dVp — F — Ft 
dt ~ ZV P' f JV dt — r P r i 

dE u _ Tp 
dt ~ "r^u 

= gV(T e -T) + qi V, 
F - -VU- 

U =U V + U em + U e + Ug, 

Ti = 7^ 7 H{Q 1 - Qi)+ 
+ q T^rH{Q l - Q 2 ) + T L H{Q % - Q x ), 
Qi = \CiVT L , Q 2 = Qi + L m -V. 

(36) 

Here (1/2)1/ is a particle volume, Qi is a thermal en- 
ergy of a particle, Tl is the melting temperature of silver, 
L m is the melting heat per Ag unit volume, H(x) is the 
Heaviside function. 

Here we assume, that the thermal capacity Ci does not 
depend on temperature. 

IV. RESULTS 

In this section only typical examples of solution of 
Equation system ([36]) will be given. A detailed analy- 
sis of these solutions for various individual parameters of 
the system, and for their arbitrary combinations will be 
presented elsewhere. 

The equations (|36p were solved numerically for vari- 
ous initial interparticle distances in bisphere, as well as 
for various intensities and frequencies of laser radiation. 
Fig. [TTb . IT2l show the calculated behaviour of the inter- 
particle distance h, ionic lattice temperatures T and elec- 
tron T e for the initial value h% = 2 nm, the laser radiation 
frequency ui = ojqi (woi is and the initial resonant fre- 
quency of the bisphere) and the pulse energy W = 10 mJ 
(the rectangular pulse). Note that the results in Fig. [TTI 
were obtained for the pulse duration r< = 10 ns (and in- 
tensity I = 10 6 W/cm 2 ) and the results in Fig. [T2l were 
obtained for n = 30 ps (I = 3.3 • 10 8 W/cm 2 ). 

As can be seen from Fig. [TTJ the particles of bisphere 
and surrounded medium are quickly heating up, and as 
a result of it the elasticity module of a polymeric adlayer 
decreases that leads to the collapse of bisphere (h — > 0) 
under the action of the Van-der-Waals forces (to ap- 
proaching of particles to minimal interparticle distance). 
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FIG. 10: Potential surface U(h,oj) (a) for a bisphere with the 
initial distance hi = 2 nm. Potential curves (b) for different 
to: 1 — 2.5 • 10 15 s" 1 ; 2 — 4.2 ■ 10 15 s _1 ; 3 — 4.45 ■ 10 15 
s -1 ; 4 — 4.53 • 10 15 s _1 . For better illustration Curves 1-4 
are shifted along the vertical axis relative each other. 



The frequency u>i (uji — u>oi(t) , ui2 = ^a2(t) for t >0) 
decreases with the interparticle distance h. Radiation 
with the frequency ioi falls out of resonance, the radia- 
tion absorption drops and the bisphere cools down. By 
the moment of collapse (t c ) the lattice temperature Tj 
becomes less than temperatures of melting T; of the par- 
ticle's material. 

That is, the resonant frequencies uj\ and <x>2 decrease 
as a result of such influence. In this case, resonant do- 
mains in a large aggregate are subjected to restructuring, 
which shows itself as photomodification of the aggregate. 
As has been mentioned before, the major contribution to 
photomodification is made by changes in the local struc- 
ture of aggregates, or, more precisely, by changes in the 
structure of resonant domains. In the approximation of 
domains by bispheres, the key role is played by the heat- 
ing of bispheres and by the Van-der-Waals forces. The ef- 
fect of light-induced forces F em = —dU em /dh on bisphere 
particles at the given radiation intensities and mean dis- 
tances (h ^ 2 nm) is negligible. Note also that at the 
given intensities the electron temperature (Fig. [Tib ) is 
just above the temperature of the lattice. 

An entirely different picture is observed in the case of 
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FIG. 11: The kinetics of h,Ti,T e under applied nanosecond 
laser pulse (/ = 10 6 W/cm 2 , r = 10~ 8 s) in a bisphere with 
initial interparticle distance hi — 2 nm (t c is the moment of 
collapse) . 



a picosecond pulse (Fig. [12"]) with n = 30 ps. As can 
be seen from Fig. [T2"r . the interparticle distance h starts 
growing as soon as the picosecond pulse is applied and it 
keeps growing after the pulse is over. This is associated 
with the light-induced multipole forces F em increasing 
with the radiation intensity and at u>i = oj\ their effect is 
to make the particles go apart (dU /dh > as shown in 
Fig. 110b . Curve 3. The particles continue moving away 
from each other for a while after the action of the pulse is 
over gradually slowing down due to friction and the Van- 
der-Waals force and eventually they start collapsing. 

Another distinctive feature of the process is a higher 
temperature of the ion subsystem which is still heated by 
electrons after the pulse due to the considerably higher 
electron temperature. As can be seen from Figs. [T2"b and 
I12[ unlike nanosecond pulses, the picosecond pulse dura- 
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FIG. 14: Temperatures and the distance h vs laser radiation 
frequency (u> — u>i) at initial value hi — 2 nm. 1 — T max , 2 
- T c , 3 — h, 4 — a{u). 



FIG. 12: The kinetics of h,Ti,T e under applied picosecond 
laser pulse (7 = 3.3 ■ 10 8 W/cm 2 , r = 3 • 10" 11 s) in a 
bisphere with initial interparticle distance hi — 2 nm. 



tion is not long enough for temperatures of electrons and 
ions to become balanced. On top of that, high tempera- 
tures at picosecond pulses are explained by the fact that 
the h distance changes very little within the pulse dura- 
tion time and uji ~ uji, which is different from nanosecond 
pulses where the bisphere collapses and falls out of reso- 
nance with the laser radiation in the middle of the pulse 
duration. 

In addition we studied the dependence of the photo- 
modification kinetics of bisphere on the frequency of laser 
radiation. The values of maximal temperature T max , as 



well as the temperature (T c ) and interparticle distance h c 
at the end of a pulse r = 10 ns and I — 10 6 W/cm 2 for 
different initial interparticle distances are shown in Fig. 
[P31 [Til We also studied dependence of the bisphere pho- 
tomodification kinetics on the laser radiation frequency. 
FiguresfT3 lll4l show the maximum temperature T max , the 
collapse temperature T c and the interparticle distances h c 
at the end of the pulse r = 10 ns and I = 10 6 W/cm 2 
calculated for different initial interparticle distances. 

As one can see from the figures, with the laser radia- 
tion frequency u>i coming close to the resonant frequency 
of the bisphere, the temperature of partricles is grow- 
ing and the distance between the particles h is reducing 
and the dependence h c (u/i) exhibits a dip. This dip how- 
ever is insignificant at large initial distances hi (Fig. 1X5]) 
and it is virtually symmetrical relative to u>qi. Such an 
insignificant change in h c is due to the Van-der-Waals 
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FIG. 15: Changes in a plasmon absorption spectrum of bi- 
sphere when particles collapse: Curves correspond to: 1 - 
the spectral band with resonant frequency woi before collapse, 
2 — the spectral band with resonant frequency u>i after col- 
lapse, Curve 3 — short- wavelengths resonant band {0J2) after 
collapse (before collapse u>2 almost coincides with w p i), Curve 
4 — the laser radiation spectrum with a frequency uii. 

force being relatively weak at large h\. The insignificant 
change of u>i does not lead to a shift of therefore heating 
of the bisphere continues throughout the pulse duration. 
The temperature maximum is achieved at the end of the 
pulse {T max =T C ). 

Quite a different situation develops in the case of small 
initial distances (Fig. fTl)) . In this case the Van-der-Waals 
force is stronger compared to the first case and the bi- 
sphere, when heated to high enough temperatures, col- 
lapses within the pulse duration. The resonant frequency 
lo\ is shifted and heating of particles ceases. Thus, the 
maximum temperature is achieved before the collapse, af- 
ter which the bisphere starts cooling down (T rnax > T c ). 

Another distinctive feature is the asymmetry of the 
temperature curves T max (uji) and T c (w>z) relative to 
T c (oji) , and also a noticeable shift of the dip in the de- 
pendence h c {u)i) . To explain this asymmetry, let us look 
at Fig. [T5l showing the bisphere absorption spectra before 
and after collapse and the laser radiation spectrum. 

For a negative initial shift Aw = loi — woi < 0, the dis- 
tance h is getting shorter under heating and the resonant 
frequency woi (^01 — * ^>i) reduces either. As a result, 
the (—Aw) shift as well as the absorption grow and the 
heating increases. That is, there is a positive feed-back 
between the radiation and the bisphere. However if the 
initial Aw, the shift of under heating leads to a lower 
absorption (a negative feed-back) and hence to a slower 
rate of collapse. 

So in the case of negative initial resonance detunings, 
Aw < 0, particles are heated to higher temperatures as 
compared to positive detunings, which accounts for the 
asymmetry of the temperature curves and the shift of the 
dip in the absorption spectrum during laser modification 



of large fractal aggregates. The change of hi as well as 
of Aw is negligible for larger initial distances h (Fig. ITS]) 
due the Van-der-Waals forces being small. This explains 
why photomodification depends so little on the sign of 
the shift. 



V. CONCLUSION 

As can be seen from the above discussion, the inter- 
action of radiation with an elementary nanoparticle ag- 
gregate (a bisphere) involves a variety of processes. A 
detailed analysis of each of them is an independent, sep- 
arate task. It is obvious that a comprehensive analysis 
of these problems is only possible on the basis of an in- 
terdisciplinary approach. Despite certain restrictions in- 
herent in our model, it enabled us to establish a number 
of fundamental features of the photomodification process 
in a bisphere as well in a fractal aggregate through the 
change of its resonant domain structure. In particular, 
the results obtained prove that photomodification of frac- 
tal aggregates due to the radiation interaction with the 
lowest-frequency resonances of the long- wavelength band 
of the plasmon absorption of an aggregate (corresponding 
to small initial interparticle gaps) is caused by shortening 
of the interparticle distance and by the associated shift 
of resonances to the long- wavelength range. In the case 
of large initial values of interparticle gaps (thick adlay- 
ers) and radiation interaction with mainly high-frequency 
resonances of the long- wavelength absorption range of a 
fractal aggregate (being closer to u> p i of a single parti- 
cle) , relative shifts of particles in a domain are insignif- 
icant even at high temperatures. Therefore the primary 
change in the domain has to be caused by the change 
in the size of particles (evaporation) down to their dis- 
appearance, and possibly by formation of plasma. As 
follows from the above reasoning, the nature of the pho- 
tomodification energy threshold both for a bisphere and 
for nanoparticle aggregates in various spectral ranges can 
be different: at small interparticle gaps it is the shift of 
particles whereas at bigger interparticle gaps it is the 
change of the sizes of particles and their disappearance. 
Another important feature established in our studies at 
small initial values of interparticle gaps is the asymme- 
try of the spectral range of bisphere collapse with respect 
to the resonant frequency of the bisphere (the radiation 
frequency causing the collapse) . As follows from Fig. [l~4l 
the bisphere collapse is mainly observed when radiation 
is absorbed in the long-wavelength range of the bisphere 
absorption curve (Curve 1 in Fig. fT5|) . If a fractal aggre- 
gate is treated as a set of resonant domains (in the ele- 
mentary case - bispheres) with various eigen frequencies, 
then collapse of bispheres will predominantly be observed 
for bispheres with frequencies exceeding w/. In that case 
a dip will appear in the inhomogeneously broadened ab- 
sorption spectrum of the fractal aggregate due to the 
laser- induced photomodification and the center of the dip 
will be somewhat shifted to the short-wavelength range 
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from loi ! 

This result of the present study is of a primary impor- 
tance, because the revealed asymmetry of the dip po- 
sition, developing wihin the laser pulse duration, can 
be one of the main factors accounting for the nega- 
tive nonlinear refraction of silver hydrosols observed in 
particular, i n 15 i 16 at the applied radiation wavelength 
A; = 1.064 /im. A detailed analysis of this result for 
nonlinear refraction and nonlinear absorption will be 
made elsewhere. However it is obvious that such spectral 
changes become one of causes responsible for nonlinear 
absorption. 

In the case of a short picosecond laser pulse, the ap- 
proach of particles turns out to be much smaller, or 
absent at all (for example, at a low intensity which is 
close to the threshold values). Thus, our analysis shows 
that under conditions similar to the experimental ones 
the optical characteristics of the resonant domain of a 
fractal aggregate are dynamically changing within laser 
pulse duration, which generates nonlincarity of the op- 
tical response of the system and is manifested in vari- 
ous nonlinear-optical processes in pulsed laser fields de- 
scribed in the Introduction. 

Our model, although a simplified one, enabled us to 
study the role of major factors and to estimate the ne- 
cessity of taking them into account. For example, with 
the particles reaching high temperatures even at mod- 
erate intensities one has to take into consideration their 
evaporation. High temperatures of electrons make it ob- 
vious that their thermoemission should also be taken into 
account. As for the polymeric adlayer properties, the 
results obtained make it necessary to allow for the ad- 
layer heterogeneity (the elasticity module of the adlayer 
increases in the direction of the particle surface). We 
could also mention a number of other problems we faced 



when using the developed model. In particular, evapora- 
tion leads to a decreased a size of particles and a reduced 
interparticle gap. The optical properties of a bisphere 
change as well as. Therefore the problem has to be con- 
sidered for a bisphere with particles of a different size. 

Besides, due to the principles incorporated in the 
model, it is possible to analyze the behavior of aggregates 
with not only metal particles but aslo with nanoparticles 
of any other material in a laser radiation field. It is quite 
obvious that the proposed model requires and what is im- 
portant allows a further improvement and development 
to provide a more adequate description of the process 
of local modification and the associated wide variety of 
nonlinear optical phenomena. We emphasize that the 
processes observed in our study have an inevitable ef- 
fect on all nonlinear-optical and photophysical phenom- 
ena which have been investigated so far but the physical 
nature of which remained obscure. One cannot rule out a 
possibility of similar processes occurring in nanocompos- 
ite materials containing metal nanoparticle aggregates in 
solid matrixes. The movement of particles in resonant do- 
mains in this case may become possible because of heat- 
ing of particles by pulsed laser radiation and melting of 
the adjacent volume of the dehydrated polymer or even 
of the adjacent glass. Other processes though associated 
with, for example, particles evaporation can also occur 
in nanocomposites, resulting in anisotropy of the shape 
of the particles. The latter as well can affect the form 
and the dynamic shift of the plasmon absorption band 65 . 
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